Local induction approximation in the theory 
of superfluid turbulence. 
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Abstract 

The local induction approximation (LIA) of the Biot-Savart law is often used for numerical and 
' analytical investigations of vortex dynamics (in particular in the theory of superfluid turbulence). 

In this paper, using renormalization group (RG) methods, some features of the LIA is considered. 



The exact statistical solution of the LIA equation is presented. The problem of "marginal" terms, 
appearing at the Wilson's approach to the RG-procedure, is concerned. 
PACS: 67.40.Vs, 47.32.Cc, 05.10.Cc 
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2 ■ 1 Introduction. 

It is accepted that chaotic vortex structures, so called vortex tangles (VT), appear in volume of superfluid 

S ■ 

^ , helium under the particular conditions. This objects are very interesting because it is well known that 

^ ■ hydrodynamic and thermodynamic properties of Hell in many respects depend on this structures [1], 

^ ' [2]. Obviously, that only statistical consideration is possible for the VT evolution. The vortex tangles 

consist of separate vortex loops, which evolve under the vortex dynamics laws (the nonlocal Biot-Savart 
, law plus interaction of vorticies with the normal component of Hell). Besides, these loops can merge and 

■ break up during temporal evolution. These two mechanisms are not independent, there is their mutual 

influence. In general, any analytical investigation of the VT evolution, based on consideration of the 
exact vortex dynamical equations is the extremely formidable task. This connected with the absence, at 
this moment, of any analytical methods for solving nonlinear, nonlocal (as the Biot-Savart), stochastical 
equations. Practically all of the known works on this theme represent various numerical (computer) 
models to explore properties of VT. In turn, the numerical calculations also have large difficults. Because 
the Biot-Savart interaction is the nonlocal one, to determ motion of a single vortex point, it is necessary 
to take into account the influence of all entire vortices to this point. Thus, it reqires unreal computer 
resources and enormous time of the calculations. 

In the initial numerical works the local induction approximation (LIA) of the Biot-Savart law was 
used [3], [4]. Some analytical advance also connected with this method [3]. Since then the vast list of 
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numerical works has been appeared. The authors obtained various resuhs, using the Biot-Savart law as 
well as the local induction approximation. However, to this day the precision and the validity of the LIA 
is the objects of discussions. 

The processes of reconnections of vortex lines even greater complicate the consideration of VT. In 
this paper we will not concern this problem, and only expose some properties of the LIA and, in general, 
RNG-based solutions of dynamical tasks. 



2 Dynamical equations. 

To describe the dynamics of a vortex filament we use the equation derived in the papers [3], [4]. To 
introduce stochastic evolution we use the Langevin approach, i.e. we add into the equation a random 
force term. That is the usual way to study turbulent- like phenomena. As the resiilt, in the local induction 
approximation the equation of motion of a quantized vortex filament in Hell takes the form; 

^^=/3Px.'' + i..'' + /(e,t) (1) 

Here t) is a radius vector of a point of a line labeled by the variable ^, s' is the derivative on 
the parameter ^; Ms time; the quantity (3 is the coefficient of nonlinearity (3 = ^log^, with the 
circulation k and the cutting parameters R (the external size, i.e. the averaged radius of curvature) and 
To (the vortex core size). The coefficient of dissipation v appears in the eq.(l) when external counterflow 
is absent [3], [4]. The additional term /(^, t) models the external random disturbances such as the "white 
noise". The concrete formulation of / will be defined later. Let's note that the eq.(l) is valid only if the 
normalization of the parameter ^ is the arclength. This is the additional approximation if to concern the 
entire consideration of vortex dynamics, since the overall length varies with time. But the aim is only to 
expose some characteristic features of the eq.(l). 

In the one- dimensional Fourier-representation the equation (1) takes the form: 

where the operator 

T120 = iP J J J J '^^i ^^^2 dwi (ici;2fciA;2e"''''^(5(fci + — fco)(5(wi + 0^2 — Wo), (3) 

and Sfc^t^, is amplitudes (in general, complex) of the Fourier harmonics of the function s{^,t), the latin 
indices denote the arguments, i.e. Sj = s{ki,cOi); the Greek superscripts means the spatial components; 
ga;/37 jg ^j^g ^j^j^ antisymmetric tensor. It is supposed, that the correlator of the random disturbances 
looks like: 

< f{kl,L0l)f{k2,L02) >= Dk-yS{kl + k2)S{iJl+LJ2) (4) 
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3 Calculations with the RG method 



Our application of the RG-procedure in principle will be similar to the usage of the one in the turbulence 
theory. The details can be found for example in [5], [6], [7]. 

Let's briefly describe the fundamental principles of the RG-method for VT. As usual, the RG- 
procedure consists of two stages. At the first stage we eliminate high (or low) harmonics from the 
consideration. Let us separate the modes Sk,u! on low and high harmonicses: s'^^ and s^^ . Here 
*fcw ~ '^fe.'^ at k £ (Ae~', A); s^^ = Sk,u> at G {K,Ae~'') and s^,s^ = in the other cases. Thus: 

Sfe,a, = S>^ + S<^ (5) 

The quantity K is the parameter of the infrared cutting. Since Sk=o = because of the closure of a vortex 
line, the value of ii' = 27r/Z; is defined by the first non-zero harmonic at the Fourier transforming of s; L 
denotes the overall length of the vortex filament. The parameter of the ultraviolet cutting A is defined 
by the vortex core size tq, and in this case A — > oo; the quantity I is some positive number, determing a 
boundary of the modes separation. 

Assume that we consider the evolution of any hamonic Sk,u> accordingly to the eq.(2) and Sk,u> is in the 
interval of low (for example) harmonics. The idea of the RG method (such as the Wilson's formulation) 
is not to consider the interaction of Sk,u> with each of harmonics in the other interval (the interval of 
high harmonics in this case). Instead of this, the RG method suggests to take into account the ensemble 
averaged influence of the whole high region. This averaging practically is executed on realizations of /. 
The interaction between harmonics in own (in this case low harmonics) interval is considered in the usual 
way. 

Concretely for a vortex loop this flrst stage of the RG-procedure is realized in the following way. We 
substitute the rel.(5) into the eq.(2) and then averaging over the ensemble of /. As the result wq obtain 
a new, modifled equation instead of the eq.(2). Note that the value A is transformed as A — > A = Ae~'. 
This new equation doesn't describe behaviour of VT on small spatial scales — the region of wave vectors 
(Ae~' -j- A) becomes unavailable. Such procedure is similar to the Kadanoff transformation [8] . If this 
new equation has the same form (perhaps approximately) as the original one and only it's parameters 
are changed, the second stage of the RG-procedure can be done. 

On the second stage we reconvert the quantity A to the original value, i.e. we are returning to the 
description on the initial scale. In the other words, we are performing the transformation A — > A = AeK 

Let's analyse some possible transformations of the parameters at the flrst stage. For example, in the 
turbulence theory the RG-procedure gives the renormalization of viscosity. This value has a positive 
addition after averaging high modes. In the other words, the influence of small spatial scales (or high 
harmonics) leads to increasing the viscosity. That, in turn, means that the dissipation of energy is 
increased. The mechanism of that phenomena is obvious. There are two processes for the evolution of 
any harmonic. The first is the thermal dissipation, the second is the nonlinear interaction with other 
harmonics, or, in the other words, the transfer of the energy of turbulent pulses across the spectrum in the 
wave- number space. (We don't take into account the external random disturbance /). Averaging high 
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harmonics lead to increasing the dissipation of low harmonics. This mean the existence of the transfer of 
the harmonic amplitudes from the low-frequency part of the spectrum to the high-frequency one. And 
just the transfer is the mechanism of the additional dissipation. Ones again emphasize: the additional 
dissipation (increasing viscosity) connected with the influence of the opposite spatial region, in the other 
words, the flux across the spectrum is responsible for increasing dissipation. And vice versa, a change 
of the viscosity under the averaging, means, in turn, the presence of the flux across the spectrum in a 
considering dynamical system. For turbulence theory this result is well known and trivial. Let's assume 
that we have solved some dynamical task and in an analogous case, we have obtained a negative addition 
to a dissipative coefficient. It must mean that the transfer across the wave- number spectrum has the 
opposite direction (i.e. from high modes to low ones). The zero addition must mean the absence of the 
transfer. 

Now let's study the eq.(l)-eq.(2) in this context. The substitution of the rel.(5) into the eq.(2), gives 
the result: 

iu:sr = r?i7 {s{<si<+ < sf<sr > + < sf>sr > + < sf>sr >> - i^k^r + /r (6) 

The brackets < ... > mean the averaging of only high harmonics. To do this averaging, turn to the 
functional integral formalism [9]. In this approach, any correlator (for example the fourth term in the 
right hand side of the eq.(6)) is represented as a product of the functional derivatives (below denoted as 
llrj "-"^ llfj^ auxiliary fields {rj , rj) from the characteristic (generating) functional of a dynamical system. 
Not going into details [9] , we write the expression to average the fourth therm: 

J dl i,,<»(l)G(») °''(1),7'5(-1)-1/2,,°(1)C(») ^^-j 

ri=fj=0 

Here G(°)(l) is the zero-order response function G(0) "^(1) = ^J^^^^, ; C(o)(l) is the zero-order correlation 
function C(0)«'5(i) = (to clarify § see the rel.(4)). 

It is well known that the Taylor-series of the first exponent gives usual Feynman's diagram series. In 
our case this takes the form: 

. /3> 7> ^_ ^ ^ r - . -O'fl^ r S S 5 

Si So J> — 



i5r]P{l)i5r]^{2) I ^ i5fi<^{Q) i5r]P{l) i5r]^{2) 

77=^=0 

The essential integrations in the last formulae are assumed. Thus, the result for < s s> represents the 
infinite series of terms, each of that is a product of C^^\ G^^^ functions (corresponding to the lines, if 
diagrams to consider) and essential integrations. Every line C^°^ or is the result of coupling the 
derivatives: {^^} or {7^77^} accordingly. Obviously only terms, containing an even number of the 
derivatives is not equal to zero. Otherwise, ry or ry remains as the efiicient in the product and the term 
equal to zero because rj = rj = 0. Besides if C(°) or G^°^ "^{a ^ (3) is presented in a product, this 
term is equal to zero, because 5"^^ is presented in the definition of G(°) "'^ and (7^°^ It is easy to 
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observe that all terms of the series for the exact correlator < s" s^^ > {a ^ (}) are identicaUy equal to 
zero. Since a ^ j3 in the external difcrcntiation, in the internal differentiation an odd number of the 
derivatives having a or /3-components are remained and thus the line G^"-* or C^"-* "'^ is necessarily 
appeared. Thus, < s" s'^ > (a ^ f3) equal to zero and this is the exact result. 

Similar argumentation leads to < s" > is equal to zero. Hence all additional terms in the eq.(6) exactly 
equal to zero. The latter means the eq.(6) doesn't change under the averaging. Thus, the transfer whithin 
the spectrum in the LI A is absent. And this result is exact. 

4 Discussion and conclusions. 

Thus we obtained that the interaction of the form: s' x s" doesn't change the spectral structure of the 
value s. In the other words, though there is the nonlinear interaction in the system, harmonics doesn't 
affect to each other. This is the extremely unexpected result. The usual scenario looks quite otherwise. 
Usually, because nonlinearity is present, interaction of only two harmonics leads to the appearencc of 
additional harmonics. They, in turn, interact between each other and with original ones. As the result, 
after some time, the full spectrum is presented in the system. In our case in the system remain only 
originally excited modes and their amplitudes arc not changed. At the same time, their phases (skc^ is 
a complex number) must be changed. Otherwise it should be full correspondense to a linear equation. 
However, it is well known the solution (that used the LIA) of the task about the decay of kink [10], that 
is not possible in a linear case. 

Obviously, it is possible to invent an infinite number of similar equations, not having a flow in the 
wave-number space. Perhaps, all of them present pure abstract interest. Nevertheless, the single practical 
note can be done. As mentioned above, the LIA is often used to approximate the full Biot-Savart law in 
numerical models of superfluid turbulence. The question about the validity of this approach was dicussed 
repeatedly and the exact opinion is absent. Let's consider this problem in the context of this work. 
Transfer within the spectrum is certanly presented in the entire Biot-Savart law, at that, probably, the 
nonlocal transfer. (Real vortex tangle is described by the Biot-Savart law). Besides, the transfer within 
the spectrum is the very important feature of nonlinear dynamical systems. For example, whether is 
it possible to imagine any approach to a turbulence theory, where the transfer of turbulent pulses from 
large scales to small ones is absent... But just the same picture takes place, when the LIA is used for 
numerical or analytical considerations of real VT. 

The last remark concerns to the usage of RG-procedures for dynamical problems. At this moment this 
question is not fully clear. Usual way emploing this method based on a perturbation theory. As a result, 
at the Wilson's formulation of RG, an infinite series of additional terms is appeared in the renormalized 
equations. So-called "marginal" terms represent, at this case, the main problem. The total value of this 
additions is uncontrolled when the Wilson's RG-procedure is carried out. The detal consideration of this 
question can be found in the paper [11]. 

In principle, we can expand s< to a perturbation series. In the other words, to employ the usual 
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approach to the RG-proccdurc. If wc have done this, the following result should be obtained. There 
is transfering amlitudes of harmonics across the spectrum is existed. And it is interesting to note, that 
the flux of harmonics is directed from small spatial scales to large ones, i.e. to the opposite side, in 
comparison with the turbulence theory. It is very unexpected result for superfluid turbulence, if it were 
right. In real, the LIA does not describe the flux of amplitudes of harmonics within the spectrum. 

Returning to the RG-method and comparing perturbation theory with the presented above exact 
solution, one can conclude that the marginal additions are really uncontrolled, even if e-expansion to use. 

The work was carried out under support of INTAS (grant N 2001-0618). 
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